arXiv: 1502.0730lv2 [astro-ph.CO] 10 Jun 2015 
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We hypothesize that cosmic string loops are the seeds about which globular clusters accrete. 
Fixing the cosmic string tension by demanding that the peak in the distribution of masses of objects 
accreting onto string loops agrees with the peak in the observed mass distribution of globular clusters 
in our Milky Way galaxy, we then compute the expected number density and mass function of 
globular clusters, and compare with observations. Our hypothesis naturally explains why globular 
clusters are the oldest and most dense objects in a galaxy, and why they are found in the halo of 
the galaxy. 


I. INTRODUCTION 

Globular clusters are spherical collections of stars 
which are found in the halos of galaxies which orbit the 
galactic center. Globular clusters are older and more 
dense than other star clusters found predominantly in 
the galactic disk. Our own Milky Way galaxy contains 
about 150 known globular clusters. Their mass distribu¬ 
tion (shown in Fig. 1) has a peak at a mass of about 
10 5 M 0 , where M 0 is the solar mass. The origin of globu¬ 
lar clusters is still not well understood (see e.g. [T] for a 
review). 

In this Letter we propose that globular clusters are 
seeded by cosmic string loops. Cosmic strings exist 
in many particle physics models beyond the Standard 
Model (see [2] for reviews on cosmic strings and their role 
in cosmology). If Nature is described by a theory which 
admits cosmic string solutions, then a network of strings 
will inevitably form in the early universe and persist until 
the present time [J. The system of strings consists of a 
network of infinite strings with a typical curvature radius 
£(t) which scales with the Hubble radius H~ 1 (t), where 
H is the Hubble expansion rate and t is physical time, 
and a distribution of string loops. At any given time, the 
distribution of string loops is characterized by a critical 
radius R c (t) which dominates the mass function of string 
loops. 

The cosmological consequences of cosmic strings are 
characterized by a single free parameters, namely the 
mass per unit length /x of the string (which is normally ex¬ 
pressed in terms of the dimensionless number G/x, where 
G is Newton’s gravitational constant). We will fix G/x 
by demanding that the mass which accretes about a loop 
with radius R c agrees with the observed peak mass in 
the globular cluster distribution. We then compute the 
predicted number density of such loops, and find that 
it agrees with the observed number density of globular 
clusters. We also compare the mass function of globu¬ 
lar clusters predicted by our model with the data and 
obtain reasonable agreement. With one free parameter 


‘Electronic address: alistair.barton@mail.mcgill.ca, 

rhb@physics.mcgill.ca, ling.lin2@mail.mcgill.ca 


we have therefore fixed the peak mass, the overall num¬ 
ber density and the mass function of globular clusters 
1 . Note that the accretion of matter about string loops 
starts as soon as loops are created, i.e. much earlier than 
the time of reionization. Hence, our mechanism offers an 
explanation for why the globular clusters are the oldest 
components of a galaxy. The string loops are initially dis¬ 
tributed throughout the region which eventually falls in 
to form the galaxy and will hence end up in the galactic 
halo rather than the disk. 

Note that in our current analysis, we are neglecting 
both the initial translational velocities of the string loops 
as well as the “rocket effect”, an anisotropic gravitational 
radiation reaction which can accelerate a string loop. 

In the following, we first review the distribution of cos¬ 
mic string loops resulting from any particle physics the¬ 
ory which gives rise to strings. In Section 3 we present 
our computations. We fix the one free parameter of the 
model, the string tension G/x, by demanding that the ob¬ 
jects accreting around string loops which dominate the 
loop mass function have the same mass as those corre¬ 
sponding to the peak in the mass distribution of glob¬ 
ular clusters in the Milky Way. We then compute the 
predicted overall number density and the mass function 
of objects accreting about such cosmic string loops and 
compare with the observed mass distribution of globular 
clusters in our Milky Way galaxy. Section 4 presents our 
conclusions. 


II. COSMIC STRING LOOP DISTRIBUTION 

Cosmic strings are one-dimensional topological defects 
which arise in a large class of particle physics models 
beyond the Standard Model. They are lines of trapped 
energy density, and this energy gives rise to interesting 
signals in cosmology. The important point [Tj is that if 


1 There is a slight caveat to this statement: the predicted number 
density also depends on constants which are known from analyt¬ 
ical considerations to be of the order 1 but whose exact values 
have to be determined by numerical cosmic string evolution sim¬ 
ulations. 
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Nature is described by a particle physics model admit¬ 
ting cosmic string solutions, then causality predicts that 
a network of strings will inevitably form in the early uni¬ 
verse and will persist to the present time. 

The energy scale 77 of the new physics which yields cos¬ 
mic strings determines the string mass per unit length /x. 
In natural units the relation is /x = ?? 2 . The amplitude of 
the gravitational effects of the strings depends linearly on 
/i. Strings come in two types: a network of infinite strings 
with a correlation length £ (t) which is proportional to the 
Hubble radius, and a distribution of string loops. Long 
string segments produce a conical discontinuity in space 
which leads to characteristic signals in CMB temperature 
maps [ 5 ], polarization maps [Bj and 21 cm redshift maps 
[7j (see [5] for a recent overview of observational signals 
of long strings). 

String loops can also play an important role in cos¬ 
mology. In fact, in early studies of the role of strings 
in cosmology it was postulated |9J that strings could ac¬ 
count for all of the observed structure. Specifically, it was 
postulated that galaxies in their entirety were seeded by 
strings. The one string - one galaxy hypothesis led to 
a value of the string tension G/x ~ ICG 6 , which corre¬ 
sponds to the energy scale 77 of particle physics “Grand 
Unification”. However, a cosmological model in which all 
fluctuations are due to strings leads to incoherent pertur¬ 
bations and predicts no acoustic oscillations in the CMB 
angular power spectrum ina- Thus, this model is ruled 
out by the CMB data 01. In fact, precision CMB ob¬ 
servations such as ACT [12] . SPT [13], WMAP [13] and 
Planck [15] lead to an upper bound on the string tension 

of m3 

G/x < 2 x 1CT 7 (1) 

(see [TJj for older results). Tighter bounds can possibly 
be obtained by analyzing CMB maps in position space 
m- Constraints on the amplitude of the stochastic back¬ 
ground of gravitational waves produced by decaying cos¬ 
mic string loops lead m to a stronger bound on G/x 
which, however, depends sensitively on the ratio of loop 
production radius to the Hubble radius (a number which 
the CMB bound is not sensitive to). 

As the discussion of the previous paragraph has shown, 
cosmic strings cannot provide the dominant source of 
fluctuations. The dominant source must be due to pro¬ 
cesses in the early universe such as inflation [ 20 ] or alter¬ 
natives (see m for a review of alternatives such as the 
“Matter Bounce” 125] or “String Gas Cosmology” [23]). 
Nevertheless, cosmic strings could have played an impor¬ 
tant role in cosmology, as long as the tension obeys the 
bound (JTJ) . For example, in many inflationary models 
cosmic strings will form at the end of inflation [23]. Fur¬ 
thermore, any phase transitions in matter which occur 
after the phase in the primordial universe which produces 
the dominant source of fluctuations (e.g. after the phase 
of inflation) will produce cosmic strings as long as the 
particle physics model admits string solutions. Hence, it 
is very interesting to search for the signals of strings in 


cosmology. If the effects are not seen to the level of a 
particular tension /x, we will have ruled out all particle 
physics models yielding strings with symmetry breaking 
scale 77 larger than ^//x. If strings are found, we can use 
them to explain cosmological phenomena which up to 
now are not explained. It is a particular application of 
the second possibility which we study in this Letter. But 
before explaining our ideas we must return to the review 
of the cosmic string scenario. 

Cosmic string loops are typically generated by the in¬ 
tersection of segments of the long strings, segments which 
move through space at relativistic speeds. Causality 
tells us that the correlation length of the string network 
(which gives the average separation and average curva¬ 
ture radius of the long strings) is bounded from above 
by the Hubble radius. Dynamical arguments show that 
the correlation length at time t cannot be parametrically 
smaller than t, otherwise string loops would be copiously 
produced and £(f) /t would increase. The network of long 
strings thus takes on a “scaling solution” according to 
which the statistical distribution of strings is independent 
of time if all lengths are scaled to the Hubble radius. 

The scaling distribution of strings implies that loops 
will be formed at any time t with a typical radius Rf(t ) 
given by 

Rf{t) = f , ( 2 ) 

where a is a constant which must be determined from 
numerical simulations. Recent simulations [251 [23: give 
a ~ 0.1. The average length of a string loop is l = /3R , 
where /3 is a constant which would be 27r if the strings 
were exactly circular. We shall use the value j3 = 10. 

Once formed, the number density of loops of radius 
R redshifts because of the expansion of space. Loops 
also slowly decay by emitting gravitational radiation m- 
Loops with radius smaller than 

R c (t) = ^G[it (3) 

will decay in less than a Hubble time. Here, 7 is an¬ 

other constant which is determined from numerical sim¬ 
ulations. Its value is of order 7 ~ 10 2 . Combining these 
results, we find that the number density n(R, t) of string 
loops of radius R is given at times t > t eq (where t eq is 
the time of matter-radiation equality) by 

n(R,t) ~ t~ 2 R~ 2 (4) 

for loops which form after t eq and 

n(R,t ) = JVa 6/2 0 -5/2 i^ 2 t -2 iT 5/2 for i? > 7 G/xt 

n(R,t ) = Afa 5 / 2 /3 -5 / 2 fg/ 2 7 -5/2 (G/x) -5 / 2 t -9 / 2 

= const for R < qG/xf , (5) 

for loops which form before t eq . In the above, the con¬ 
stant N depends on the average number N of long string 
segments per Hubbe volume. Since two long string seg¬ 
ments are required to form a loop, we could expect N 
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to depend on the square of N. We will be interested 
in values of G/x for which the critical radius © at the 
present time to is smaller than aP~ 1 t eq . Such loops form 
before the time of equal matter and radiation, and their 
distribution is given by ©>• 

The total number density of strings is obtained by in¬ 
tegrating <[5j) over the radius. The integral is dominated 
by R ~ R c (t). For times between t eq and the present 
time to we obtain 

ntotalit) ~ 7Va 5 / 2 /3- 5/2 7 - 3/2 2- 3/4 (G/x)- 3/2 r 3 , (6) 

where z(t) is the cosmological redshift at time t, and 
where the constant 5/3 resulting from integrating over 
the loop distribution ([5| has been absorbed in the factor 
N. Note, in particular, that the number density is larger 
for smaller values of G/x. Hence, for such values of G/x 
there will be many string loops embedded in the region 
which collapses to form a galaxy (the galaxy seed being 
given not by cosmic strings, but by the dominant source 
of fluctuations). In fact, the accretion of matter about 
the galactic center sweeps up string loops from a larger 
region of space. Working in the Zel’dovich approximation 
[25] , the local number density of string loops inside a 
galaxy will be enhanced by a factor of F. Accretion and 
virialization in each direction leads to the estimate F = 
64. 

In the following we will explore the possibility that 
the string loops are the seeds for the globular clusters 
embedded in the galactic halo. 


III. GLOBULAR CLUSTERS FROM COSMIC 
STRING LOOPS 

According to linear cosmological perturbation theory, 
accretion of matter about a cosmic string loop starts 
(modulo logarithmic growth at earlier times) at t eq . At 
that time, the mass function of string loops is dominated 
by loops with radius 

R = ~(GlJit eq . (7) 

Whereas these loops will have decayed by the present 
time, the objects they seed will keep growing. The local 
number density of such objects inside a galaxy is hence 
given by 

= FNa^r^tlft^Rcite,)- 312 

= FNa 5/2 p- 5/2 j- 3/2 zV 2 {GtJ,)- 3/2 to 3 . ( 8 ) 


The mass which has accreted about these seed loops at 
the present time 2 is given by 

M(jGnt eq ,to) = 

= £ 7 (G/x) 2 0 - 1/2 |. (9) 

Inserting the values of to and G we find 

| ~ 10 23 M o . ( 10 ) 

Making use of the values /3 = 10, 7 = 10 2 and z eq = 10 4 
and setting the mass © equal to the peak of the mass 
function of globular clusters in our galaxy Mqc ~ 10 5 M 0 
we find 

G/x ~ 10 -9 ' 5 . (11) 

Inserting this value into the local number density (| 8 |, 
using (in addition to the previously mentioned values for 
p and 7 ) the values of the constants (see [26]) a = 0.3, 
N = 10 2 , and inserting the value of to which yields 3 ~ 
10 ” 21 kpc -3 we find 

niocai ~ 10 ~ 2 (kpc)~ 3 . (12) 

This is of the same order of magnitude as the observed 
number density of globular clusters inside the Milky Way 
galaxy. 

The cosmic string loop seed hypothesis for global clus¬ 
ters predicts the primordial mass function of globular 
clusters. The value of M c , the peak in the mass func¬ 
tion, was used in our analysis to fix the string tension. 
However, once the tension is fixed there is no more free¬ 
dom in the mass distribution. For M > M c the mass 
function scales as M -5 / 2 (as follows directly from the 
string loop distribution ([ 5 ])). For masses smaller than 
M c we predict a linear decay. This follows from the fact 
that the loop radius distribution is flat but loops with 
radius smaller than yG/xt live only a fraction of a Hubble 
time step which scales linearly with R. 

In Fig. 1 we compare the predicted mass function (the 
solid line) with the observed distribution of globular clus¬ 
ters in the Milky Way (histogram values). The values for 
the blue histogram boxes come from the compilation in 
[3] of the properties of globular clusters in the Milky Way. 
The theoretical curve is obtained by taking the comov¬ 
ing number density n(R , t eq )z 3 q of loops (where n(R , t) is 
given in ©). as a function of R , multiplying the result by 
the concentration factor F , allowing each loop to grow in 
mass by a factor of z eq (independent spherical accretion), 
and converting the distribution into a mass distribution 
n(M,to), taking into account the Jacobian of the trans¬ 
formation from R to M. The result is then multiplied 
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Wlocal 


We are here assuming that the accretion continues to the present 
time. 
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FIG. 1: Mass function of globular clusters in the Milky Way galaxy. The horizontal axis is mass on a logarithmic scale, the 
vertical axis gives the number density on a linear scale. The histogram show the data taken from [3]. The solid curve is the 
prediction of our model for the value of G/z = 5.43 x 10 -10 which minimizes y 2 . The cosmic string parameters chosen are 
described in the text. Note, in particular, that the overall amplitude of the theoretical curve is predicted once G/z has been 
fixed to fit the peak location. 


by the bin size SM and by the volume of the Milky Way 
galaxy. Taking the bin size to be SM = fM , where / is 
a number, we obtain for the peak number density bin 

SN = NFfa 5/2 p- 5/2 'y- 3/2 zl / q 2 

x(G/r)- 3 / 2 yto 3 (13) 

where V is the volume of the galaxy. The peak mass of 
the theoretical curve is given by In Fig. 1 the hori¬ 
zontal axis is mass (logarithmic scale) while the vertical 
axis gives the number of globular clusters (linear scale) 
inside of the Milky Way galaxy. The exact value of G/z 
for the theoretical prediction has been tuned to give the 
lowest value of y 2 between the theory curve and the data. 


IV. CONCLUSIONS AND DISCUSSION 

We have seen that both the number density and mass 
of string loops with a mass per unit length given by © 


Its value is G/z = 5.43 x 10 -10 . 

In Fig. 2 we show how the peak position and amplitude 
of the mass function shift as G/z varies. As G/z decreases, 
the peak mass decreases but the predicted number of 
globular clusters increases. The surprising result is that 
a good fit of both amplitude and peak number density 
occur for the same value of G/z 3 . 

Our analysis does not take into account interactions 
between globular clusters and other processes which can 
change the mass of a globular cluster as a function of time 
during the highly nonlinear process of galaxy formation. 
Another effect which we have neglected in the current 
analysis is that of initial translational motion of string 
loops. 


is correct to explain the origin of globular clusters. We 
hypothesize that globular clusters are seeded by string 
loops which dominate the loop mass distribution at the 
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FIG. 2: Dependence of the mass function of our model on the string tension Gfj,. The curves shown are for Gp = 2.92 x 
10~ lo ,G/x = 3.98 x lCP 10 ,G/x = 5.43 x 10~ 10 and G/x = 7.41 x 10~ ll] (in increasing order of mass at the peak position). The 
axes and data are the same as in the previous figure. 


time t eq when the accretion of cold matter onto the loops 
can start. We emphasize that fixing one free parameter 
has allowed us to fit both the mass and the number den¬ 
sity of globular clusters. In addition, the model predicts a 
mass function which is in reasonable agreement with ob¬ 
servations. Our model for the origin of globular clusters 
also naturally explains the fact that the globular clus¬ 
ters are the oldest star clusters in a galaxy, that they are 
the densest regions, and that they are distributed in the 
galactic halo rather than confined to the galactic disk. 

Our theory predicts the the number of globular clusters 
in a galaxy will scale linearly with the mass of the galaxy. 
It also predicts that the mass function of globular clus¬ 
ters should be universal across galaxies. Since the energy 
distribution of dark matter falling onto a cosmic string 
is more peaked than that falling into a Gaussian fluctua¬ 
tion [29] , we predict falling velocity rotation curves about 
the center of a globular cluster. Not all string loops are 
swept up into galaxies. Hence, our model also predicts 
the presence of compact objects of globular cluster mass 
in the field between galaxies. The corresponding number 
density does not include the factor F ~ 10 2 . 

In the current analysis we have neglected the initial 


translational motion of the string loops. Since this pecu¬ 
liar velocity redshifts and since we are considering string 
loops produced long before the time of equal matter and 
radiation, the loop motion should be a small effect by the 
time that the loop starts accreting matter at t eq . We are 
currently in the process of working out the details. An 
effect which might be more important (and which we are 
neglecting here) is the rocket effect. 

A question which we have not addressed here is the 
issue of star formation inside the object which accretes 
onto the string loops (see e.g. UBECEI for recent studies of 
this question in the context of early structure formation 
from strings). There may well be a significant difference 
in the star formation processes for matter accreting onto 
a string loop inside a galaxy versus in the field. 
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